Vertex corrections to the mean-field electrical conductivity in strongly disordered 

electron systems 
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Mean-field theory ol non-interacting disordered electron systems is widely and successfully used to 
describe equilibrium properties of alloys in the whole range of disorder strengths. It, however, fails 
to take into account effects of quantum coherence and localizing back-scattering effects when applied 
to transport phenomena. We present an approximate scheme extending the mean-field theory for 
one-electron properties in that it offers a formula for the two-particle vertex and the electrical 
conductivity non-perturbatively including the leading-order vertex corrections in a way that the 
approximation remains consistent and the conductivity non-negative in all disorder regimes. 
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INTRODUCTION 



Fluctuations of the spatial distribution of the atomic 
potential give origin to the zero-temperature resistivity of 
crystalline solids. In order to study in detail the impact 
of fluctuations of the atomic potential on charge trans- 
port it is appropriate to neglect all less important agents 
influencing the response of the electron gas on external 
electric perturbations and to investigate only scatterings 
of free electrons on random impurities. It is usual to 
study the effects of impurity scatterings on an Anderson 
model of disordered electrons^ with a site-independent 
distribution of the atomic potential. A straightforward 
way to investigate this model is to use a perturbation (di- 
agrammatic) expansion in the random atomic potential. 
Nontrivial and physically interesting results can be, how- 
ever, reached only via non-perturbative approaches. The 
most reliable non-perturbative method of describing the 
effects of disorder on one-electron functions is a mean- 
field, coherent-potential approximation (CPA)<2 The co- 
herent potential approximation is nowadays considered 
as the archetype of mean-field theories of quantum disor- 
dered and interacting systems. Its generalized forrn^ of- 
fers one possible interpretation of equations of motion in 
the dynamical mean-field theory (DMFT). 4 The CPA has 
proved reliable to produce an accurate equilibrium elec- 
tronic structure of disordered systems^ as well as trans- 
port properties of random alloys in a wide range of disor- 
der strength^ It, however, fails to account for inter-site 
quantum coherence and backscattering effects. The CPA 
is essentially unable to go beyond the semi-classical de- 
scription of transport properties qualitatively captured 
by the Boltzmann equation. This inability is due to the 
fact that the CPA does not include vertex corrections to 
the one-electron (Drude-type) electrical conductivity! 7 - In 
the simplest single-band tight-binding model the electri- 
cal conductivity at zero temperature in a direction can 
be represented in the CPA via the averaged one-electron 



k, z + 



k + q, 2_ 



k' + q, z_ 



FIG. 1: Attribution of complex frequencies z\,Z2 and mo- 
menta k, k',q in the electron-hole diagrammatic representa- 
tion of the generic two-particle functions X kk / (q; zi, z-i) 
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where v a (k) = de(k)/dk a is the electron group velocity 
in direction a and superscript R denotes the retarded 
part of the averaged one-electron propagator. All the 
contributions to the electrical conductivity beyond cr' ' 
originate from coherent propagation of an electron and a 
hole and are contained in vertex corrections. 

Vertex corrections to the Drude conductivity from 
Eq. ([T]) are important in various situations. In low 
dimensions (d < 2) or for sufficiently strong disorder 
they are responsible for vanishing of diffusion called An- 
derson localization. ? Further on, tunnel conductance or 
transport through multilayered strongly disordered al- 
loys and dirty metals are essentially influenced by vertex 
corrections i^ilS. To obtain more realistic results for the 
electronic transport in dirty metals one has to go beyond 
the standard mean-field conductivity and to develop ap- 
proximations containing spatial quantum coherence and 
backscattering contributions. 

To calculate the full electrical conductivity from the 
Kubo formula with the current-current correlation func- 
tion we need to use the averaged two-particle Green func- 
tion at the real energy axisAi It can be represented as 

G&,(q,w) = G£i,{q;E + uj + iaO + ,E + ibQ + ) (2) 

with a double superscript ab specifying complex half- 
planes in which the limit to the real axis of the two energy 
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(frequency) variables is taken. We use a, b = ±1 or equiv- 
alently a, b = R,A for the upper, lower half-plane for 
the electron and hole, respectively. Fermionic variables 
k, k' are the incoming and the outgoing particle momenta 
and q is the difference between the incoming momenta 
of the particle and the hole in the pair, confer Fig. [TJ 
Since the frequencies are not dynamical variables for dis- 
ordered noninteracting systems we often suppress them 
in order to simplify notation. Actually, the unspecified 
frequency (energy) variables are set to the Fermi energy 
as all calculations are done at zero temperature. 

The contribution to the electrical conductivity from 
the free uncorrelated propagation is the Drudc term in 
Eq. (JTJ) . The contribution from the vertex corrections 
is expressed via the two-particle vertex defined from the 
averaged two-electron Green function as 



of the vertex corrections to the electrical conductivity 
derived from the asymptotic limit to high spatial dimen- 
sions in Refs. flMTi in that we keep the full structure of 
the Brillouin zones of low-dimensional lattices and hence 
treat momentum convolutions there more realistically. 



II. ELECTRICAL CONDUCTIVITY FROM A 
SYMMETRIC BETHE-SALPETER EQUATION 

To derive a non-perturbative expression for the elec- 
trical conductivity we use a representation of the full 
two-particle vertex via a Bethe-Salpeter equation in the 
electron-hole channel. Using the allocation of variables 
of two-particle quantities from Fig. Q] and denoting the 
irreducible vertex in the electron-hole channel A we have 



^kk' (l) — ^k^k+q 

x [5(k-kO+r^(q)G£,Gt +q ] • (3) 

The vertex corrections to the zero-temperature conduc- 
tivity then readi^ 
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We suppressed the bosonic momentum variable q = 
and frequency lu = as we will do throughout the paper, 
when no confusion may arise, for variables the values of 
which are set to zero. 

Troubles with the decoupling of the full conductiv- 
ity into the Drude term and vertex corrections emerge 
in the limit of strong disorder and low dimensions where 
the two-particle corrections with negative sign may over- 
weigh the positive one-electron contribution. There is not 
a universal guarantee that the total conductivity stays 
non-negative in the whole range of the disorder strength. 
Even systematic corrections to the mean-field conductiv- 
ity derived from the limit to high spatial dimensions may 
lead in three-dimensional models to negative values of 
conductivity 13 ! 14 Since the mean-field approximation is 
consistent for one-electron properties in the whole range 
of the disorder strength, it is desirable to have an anal- 
ogous consistency also for transport properties and elec- 
trical conductivity with leading-order vertex corrections 
applicable on realistic random systems in three dimen- 
sion. 

It is the aim of this paper to present a scheme how to 
reach an approximate expression for the electrical con- 
ductivity in disordered noninteracting systems extending 
the mean-field expression by involving vertex corrections 
so that to keep the conductivity non- negative. We use the 
asymptotic limit to high spatial dimensions to derive an 
approximate formula for the electrical conductivity with 
vertex corrections delivering non- negative results. The 
present paper is an improvement upon the calculation 
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This Bethe-Salpeter equation can be extended to the full 
two-particle Green function by using the definition in 
Eq. ([2]). The two-particle Green function can then be 
determined from the electron-hole irreducible vertex as 



^kk;' (q) — G^G^ 



^ k - k ') + A ?E<''( c i) G £''k'(q) 



(6) 



This is an integral equation in momentum space with a 
kernel A£ k , (q)G£, G k / +q where fermionic momenta k and 
k' are the active ones. This integral kernel is asymmet- 
ric in their active variables even for the difference mo- 
mentum q — 0. We can, however, transform the above 
Bethe-Salpeter equation to a symmetric one allowing for 
an easier diagonalization and non-perturbative solution. 
To do so we redefine two-particle functions in momen- 
tum space by multiplying them from right and left by 
one-electron propagators 
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The one-electron propagators are now absorbed in the 
rescaled vertex that forms a symmetric kernel for a 
Bethe-Salpeter equation for the two-particle vertex 

f k^ (q) = a& (q) + ^ E A k b k" (q)P^'k' (q) • (8) 

k" 

Analogously to this Bethe-Salpeter equation we obtain 
a symmetrized Bethe-Salpeter equation for the full two- 
particle Green function the solution of which we formally 
represent as 
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where * stands for the appropriate matrix multiplication 
in the eh scattering channel used in Eq. ©. 

This representation of the two-particle Green function 
can be used to derive a non-perturbative representation 
of the electrical conductivity that is free of the decom- 
position into the Drude term and vertex corrections. We 
obtain for the static optical conductivity at zero temper- 
ature 
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The expression for the electrical conductivity via the 
electron-hole irreducible vertex A in Eq. (fT0|) leads to 
non-negative results if the inverse operator on its right- 
hand side is correctly calculated or approximated. Its 
perturbation expansion, or separating the even part of 
the inverse operator, leads to decomposition of the con- 
ductivity into the Drude term (even contribution) and 
the vertex corrections (odd contribution). It is hence 
mandatory to use a non-perturbative evaluation or ap- 
proximation for the inverse operator in Eq. (|10[) to guar- 
antee non-negative results for all strengths of disorder. 



III. EXPANSION AROUND MEAN FIELD 

The representation from Eq. (|10[) offers a way how to 
calculate the electrical conductivity non-perturbatively 
without the necessity to decompose it into the mean- 
field one-electron contribution and the two-particle ver- 
tex corrections. To guarantee non- negative results for the 
conductivity one has to resolve the inverse in momentum 
space non-perturbatively. Ideal would be to diagonalize 
the irreducible vertex exactly. It is, however, possible 
only numerically in a discretized space with only a few 
points in the Brillouin zone. In realistic situations one 
need not evaluate the inverse operator exactly, since the 
vertex corrections are in most situations small. So that 
to keep the effort in calculating the electrical conductiv- 
ity with vertex corrections within reasonable limits, we 
resort to a non-perturbative approximation of the matrix 
inversion in Eq. (|10p . We use the asymptotic limit to high 
spatial dimensions as a guide in the selection of the ap- 
propriate manageable approximate scheme. The vertex 
corrections can be systematically controlled in this way. 



A. Leading-order non-local electron-hole vertex 

To determine the leading-order vertex correction to 
the electrical conductivity we use an expansion around 
the mean- held solution, being exact in d = oo. The 
mean-field solution becomes exact when the off-diagonal 
part of the full one-electron Green function vanishes. It 



leads to a local self-energy S(z). A small parameter for 
the expansion around the mean-field solution is the off- 
diagonal mean-field Green function that can be repre- 
sented in momentum space as 
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with £ = z — S(z). The mean-field approximation is 
consistent only when the one-electron Green function is 
everywhere replaced by its local element. It applies in 
particular on two-particle quantities. The only consistent 
and unambiguous two-particle vertex within the mean- 
field approximation is the local one 
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where we used the two-particle irreducible vertex A. It 
is related to the one-electron self-energy via a Ward 
identity.—^ The Ward identity is channel-dependent. We 
have in the electron-hole and electron-electron channels^ 
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%G R ~ x RA (0) 
E w z r 

qw ~ x RR(o) 
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where we denoted £ 

dG R {uj)/duj\ u=0 , and Z R = E R '/(E W - 1). We also re- 
lated the irreducible vertices with the homogeneous part 
of the non-local two-particle bubble, a convolution of two 
one-electron propagators 
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(C, C'; q) = n- 1 ]T G(k, C)G(k + q, CO 



(14) 



This bubble is non-local even in the mean-field limit (d = 
oo). 

When trying to extend the mean-field approximation 
also to the non-local two-particle functions, the results 
are no longer consistent, since the non-local vertex de- 
rived with the aid of the Ward identity, used in the 
CPA, does not cover all the non-vanishing contributions 
to two-particle functions in d — oo. This ambiguity is a 
consequence of the fact that two (linear) operations do 
not commute, namely (functional) derivative of Green 
functions with respect to external sources and the limit 
d — > ooi^Z If we introduce the off-diagonal (non-local) 
part of the two-particle bubble 
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we can represent the leading non-local vertex correction 
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r kk <(C+,C-;q) 



= 7(C+,C-) + A(C+,C-) 

7 (C + ,C-)x(C+,C-;q) 



.i-A(C+,C-)x(C+,C-;q) 
7(C+,C-)x(C+,C-;q + k + k') 



l-A(C + ,C-MC +> C-;q + k + k') 



(16) 



Standardly the full momentum-dependent mean-field 
vertex contains only the first two terms on the right-hand 
side of Eq. (|16p independent of the fermionic momenta 
k and k'. Due to symmetry reasons they do not con- 
tribute to the electrical conductivity. It is the third term 
that generates the leading-order vertex corrections to the 
electrical conductivity and is responsible for the so-called 
weak localization^ 

To avoid addition of the vertex corrections to the 
mean-field conductivity we use the Bethe-Salpeter equa- 
tion ^ and express the conductivity via the electron- 
hole irreducible vertex A. Since we expand around 
the mean-field solution, we must use a modified Bethe- 
Salpeter equation suppressing multiple scatterings on the 
same site so that to avoid double counting of scattering 
effects. We denote the modified irreducible vertex A. 
A modified Bethe-Salpeter equation for the full vertex 
read$i2 

Ikk' (z+ , z- ; q) = A kk . (z+ , z_ ; q) 

^ £ Akk<< (*+, z- ; q) G + (k")G_ (k" + q) 
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(z+,z-;q) . (17) 



The leading order contribution to the off-diagonal part of 
the irreducible vertex in high spatial dimensions is just 
the last term on the right-hand side of Eq. (jTHJ) . We then 
have 
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(18) 



To calculate the conductivity with the aid of vertex A 
and a symmetrized Bethe-Salpeter equation we redefine 
the two-electron Green function 



Ckk' (q) — G k 



+q 



5(k-k') + G:rf k ,(q)G k , +q G k , (19) 
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that does not change the long-range diffusive behav- 
ior of the full two-particle Green function. This con- 
strained function can be represented as a solution of a 
symmetrized Bethe-Salpeter equation 



Gkk'(q) = G k+q 1-A afc (q)* 
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(20) 



with a kernel that we denote for further calculations as 
K' 1 = 1-2. 



We use the constrained two-electron Green function 
G in the calculation of the electrical conductivity. The 
result is not changed by this substitution. From the sym- 
metrized Bethe-Salpeter equation (fT9"|) we obtain 
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The off-diagonal propagator G is the fundamental pa- 
rameter in the expansion around the mean-field limit. It 
makes this expansion consistent and moreover it makes 
the calculation of corrections to the mean-field result nu- 
merically more stable. It is preferable to use the full local 
mean-field vertex r y ab instead the irreducible one A afc in 
all formulas of the expansion around mean field, since 
the latter contains a pole in the RR {AA) channel that 
is cancelled in the former one. Note that the leading- 
order vertex corrections calculated from an expansion of 
the right-hand side of Eq. (|2ip coincide with the lead- 
ing corrections to the mean-field conductivity derived in 
Ref . M- 

The explicit representation for the integral kernel 
from Eq. ([2"Tj) to be inversed is 



(K a %l (q) = (L a %l (q) - 7 afc G a (k)G"(k' + q) 



'kk' 

= <5k,k< - 7 Qb G a (k)G fc (k' + q)S ab (k + k' + q) 

- 7 ab G a (k)G fc (k' + q) (22) 
with k and k' as active variables. We denoted 

s ab (Q)=r b T b (Q)/[i-j ab T b (Q)] ■ (23) 

One can try to invert numerically exactly the matrix 
from Eq. (|22[) . but it is a rather demanding task. Since we 
expect that the mean-field conductivity dominates within 
the energy bands with relevant vertex corrections only 
close to the band edges we can resort to an approximate 
matrix inversion so that non-negativity of the result is 
preserved and the bulk conductivity not much affected. 

B. Approximate solution 

The dominant contribution to the irreduciblle vertex 
A from Eq. (|18|) comes from the pole in S RA (Q) for Q = 
0. It means that we can single out fermionic momenta 
k, k' in the inversion so that k+k' + q s» 0. Then the pre- 
factor G k G k < +q at function S RA (k + k' + q) in Eq. I[2"2"|) 
becomes relevant only in the vicinity of the pole, that is 
for k + k' + q ps 0. We further get rid of the dependence 
of the pre-factor on the fermionic variable in that we 
replace it by its average over the Brillouin zone, namely 
A^ 1 G k G_ k . This simplification leads to the inverse 
of a reduced matrix 



L kk'(q) = <*kk' 



7x(0)S(k + k' + q) 



(24) 
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Although this approximation is justified for vertex A , 
since only the electron-hole channel contains the pole, we 

i • i r -r RR 

use it, due to consistency, also for vertex A 

Inverse of the matrix on the right-hand side of Eq. (|24[) 
can be explicitly evaluated by using a Fourier transfor- 
mation diagonalizing convolutions (correlations) result- 
ing from multiplication in the matrix inversion in mo- 
mentum space. We introduce a Fourier transform in a 
d-dimensional momentum space as follows 



/(x) = fdq e !q x /(q) 



/(q) = 



dx 

(2tt) c 



e" lqx /(x) 



We can now explicitly represent the inverse matrix 
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e -i(k— k')-x _|_ ^ab—ab gab ^_^ e -i(k+k' 



i - (T b T ) S ab (x)S ab (-x) 



(25) 



as well as 



conductivity in a form guaranteeing its non-negativity 
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{{vG A L RA vG R )-^{vG R L RR vG R )} . (28) 



Note that vG{x) = -vG(-x) while S(x) = 5(-x). It is 
the first term on the right-hand side of Eq. (|28p that is 
dominant within the band with non-zero imaginary part 
of the self-energy and a RA > \a RR \ > 0. 

The above approximation is consistent if 

G A L RA G R ^ (q) remains positive for q = 0. It is 

fulfilled if a stability condition 
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is satisfied. This poses a restriction on the disorder 
strength measured by the full local mean-field vertex 7 
for which this approximation is consistent and realiable. 
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where we denoted 
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We suppressed the frequency variables, the (infinitesi- 
mal) imaginary parts of which are indicated by the su- 
perscripts. 

We further introduce Fourier transforms of velocities 
as vG (x) = J dkw(k)G^e lk x so that to reach a rep- 
resentation of the electrical conductivity containing the 
leading-order vertex corrections to the mean-field Drude 



IV. ELECTRON-HOLE CORRELATION 
FUNCTION AND DIFFUSION POLE 

The proposed approximate scheme of inverting the 
kernel in the Bethe-Salpeter equation so that the elec- 
trical conductivity never goes negative can also be used 
to determine the electron-hole correlation function. Since 
we are interested only in the diffusion pole of this function 
we redefine it here so that to keep the polar structure but 
to make the explicit expression in the expansion around 
the mean-field limit as simple as possible. We hence use 
the off-diagonal one-electron propagators to define a con- 
strained electron-hole correlation function 
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As is well known this function contains the so-called diffu- 
sion pole when the Ward identity between the self-energy 
and the irreducible vertex is obeyed for all frequencies.— 
Its mean-field, local version is expressed in Eq. (fT3| . The 
full, frequency-dependent Ward identity holds only in the 
mean-field limit. Once we go beyond and introduce non- 
local corrections to the irreducible vertex, the full Ward 
identity is in conflict with analyticity of the self-energy 
and integrability of the diffusion polei 21 i 22 When we ap- 
proximate the non-local part of the irreducible vertex 
from which we then determine the self-energy, the Ward 
identity can be obeyed only in the static limit, that is, 
for zero difference frequency between the electron and the 
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hole. The exact low-energy asymptotics of the electron- 
hole correlation function with the static Ward identity 
reads^S 



^RA 



-iA E uj + D° e (oj)(i 2 



(31) 



where Ae > 1 is a constant directly proportional to the 
disorder strength. This electron-hole correlation function 
contains the diffusion pole of order for q = 0, but its 
weight decreases with increasing disorder strength and 
vanishes in the localized phase.— Physical consistency 
dictates that this low-energy behavior of the electron- 
hole correlation function should be reproduced also in 
the proposed approximation. The electron-hole correla- 
tion function in this approximation has an explicit rep- 
resentation 



-RA 



(q) = 



(q) 



I ~ 1 RA {G A L RA G R 



(q) 



(32) 



In an approximation with a non-local irreducible ver- 
tex we must go beyond the local mean-field self-energy 
to guarantee the low-energy asymptotics of Eq. ([3Tj) . 
To make the proposed approximation consistent with 
the diffusion pole we have to correct the one-electron 
self-energy appropriately. We can use the static (zero- 
frequency) Ward identity, that complies with this type 
of approximation, to determine the imaginary part of the 
self-energy from the irreducible vertex 
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^A RA ,%G R 
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where both the two-particle irreducible vertex and the 
averaged one-electron propagator are unrestricted. The 
real part of the self-energy is then calculated from the 
Kramers-Kronig relation.— 

This way of making the approximation with a non- 
local irreducible vertex compatible with the diffusion pole 
in the electron-hole correlation function is rather tedious 
and numerically demanding. Non-local parts in momenta 
and frequency of the self-energy are interconnected and 
must be determined simultaneously. There is, however, a 
more effective way how to assure a qualitatively correct 

low-energy asymptotics of function $ from Eq. (|32|) . 
Instead of correcting the self-energy we can rescale the 
local static vertex by a positive constant so that the dif- 

fusion pole in function $ E (q) is guaranteed with the 
local mean-field self-energy. For this purpose we intro- 
duce a positive number < <fi < 1 with which we rescale 
the static part of the mean-field vertex 
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in Eq. (|22p . Note that vertex 7 is not rescaled in the 
definition of function S(q) in Eq. (|23|) . since it is consis- 
tently defined within the mean-field approximation. The 



scaling parameter is determined so that the denominator 
vanishes in the static and homogeneous limit. That is 
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It is easy to rewrite this equation to a more explicit 



one 
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the solution of which can be straightforwardly found nu- 
merically. The introduced static scaling factor effectively 
decreases the disorder strength and prevents the system 
from undergoing the Anderson localization transition to 
the phase without macroscopic diffusion. This simpli- 
fied approach to make approximations free of unphysical 
behavior at any disorder strength is not good enough to 
take into account spatial coherence of electron scatterings 
on random impurities beyond the weak-localization effect 
covered by the leading-order vertex corrections. To reach 
the Anderson localization transition one has to introduce 
a self-consistent renormalization of the irreducible vertex 
Aj2£ Nevertheless, the proposed approximation with non- 
local vertex corrections to the two-particle irreducible 
vertex does not result in unphysical behavior. 

With the diffusion pole in the electron-hole correla- 
tion function we can determine its low-energy behavior 
expressed in constants Ae and D E . They are 
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1 ra {E + u,E)(g A {E) 



xL RA (E + uj, E)G'"{E + uj) 



(37a) 



-^ HA (E,E) 

xV? (G A (E)L RA (E,E)G R (E))(q 



q=0 



(37b) 



where we explicitly specified the energy variables in the 
one and two-electron functions so that to avoid confusion. 
In the mean-field solution with the local irreducible ver- 
tex where the Ward identity is obeyed for all frequencies 
the diffusion constant is related to the static optical con- 
ductivity at zero temperature as a — e 2 ppD Fl where pF 
is the density of states at the Fermi energy^ Since the 
theories with non-local irreducible vertex cannot comply 
with the Ward identity beyond the static limit, the dif- 
fusion constant and the conductivity will no longer obey 
this relation. 
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FIG. 2: (Color online) Density of states of a symmetric binary 
alloy (c = 1/2) on a simple-cubic lattice for various strengths 
of disorder. Split band occurs at A = ^Jw where w = 6 is the 
band width. 



V. RESULTS 



FIG. 3: (Color online) Conductivities of a symmetric binary 
alloy on a simple-cubic lattice calculated from the Drude for- 
mula, a' ', Drude with leading-order vertex corrections, ctq, 
and conductivity a a from Eq. ([28} as a function of disorder 
strength A. 



We use the developed approximate theory to evaluate 
the impact of vertex corrections on the electrical con- 
ductivity of a binary alloy on a simple-cubic lattice with 
a dispersion relation e(k) = — (cosfc^ + cosk y + cosfc z ). 
The site-independent distribution of the random poten- 
tial is p(V) = (1 - c)5{V - A/2) + cS(V + A/2). We use 
the mean-field approximation (CPA) on the self-energy 
and the one-electron Green function. We use elliptic in- 
tegrals to calculate integrals with the density of states 
of a simple-cubic latticed Dependence of the density 
of states on the disorder strength is plotted in Fig. [2j 
The small imaginary part of energy to regularize inte- 
grals along the real axis was chosen ij = 10 -8 . 

We need direct momentum integration to calculate 
two-particle properties in this approximation. We use the 
standard technique^ to discretize the irreducible Bril- 
louin zone of the single-cubic lattice with 8 3 to 24 3 mesh 
points according to the precision needed. All calculations 
were performed at zero temperature for a simple-cubic 
lattice and a symmetric (c = 1/2) binary alloy with a 
split-band limit for A c = \/6 at which the imaginary 
part of the self-energy diverges. 

We first compare the results for the electrical con- 
ductivity calculated from different mean-field approaches 
with vertex corrections. In particular, we compare 
the approximate zero-temperature conductivity from 
Eq. ((11} with that from the Drude formula, Eq. Q 
and with the conductivity including vertex corrections 
from Eq. (jj) using the vertex function from the high- 
dimensional asymptotics, Eq. (JTBJ) . We plotted in Fig. [3] 
behavior of these conductivities near the internal band 
edge (split band). We set the electron charge e = 1. 
We can see that before the split band limit is reached, 
the effect of disorder becomes so strong that the vertex 



corrections, when added linearly, over weigh the mean- 
field result and turn the full conductivity negative and 
the approximation inapplicable. The approximate non- 
perturbative inclusion of vertex corrections from Eq. (|28[) 
remains positive up to the band edge. Deeper within the 
band the two approximate formulas coincide, since the 
relative weight of the vertex corrections is small. 

The ratio of the modulus of the vertex correction 
and the Drude conductivity in shown Fig. |4] We com- 
pared the leading-order vertex correction derived from 
conductivity a a calculated from formula (|2"8)l . that is 
Acta = (ta — o-(°\ We also added the leading-order vertex 
corrections in high spatial dimensions Acr,^ from Ref. [TH 
and omf a mean-field conductivity with the vertex cor- 
rections from Ref. [13 where integration over momenta 
was simplified via the limit to high spatial dimensions. 
The approximation from Ref. [T3 is similar to the non- 
perturbative approach of this paper and was derived for 
the same purpose: to produce non-negative conductivity. 
It completely eliminates, unlike the present one, the di- 
vergence in the electron-hole vertex function due to the 
simplification of momentum convolutions used there. We 
can see that the latter vertex correction shows a different 
trend from the other two, that is decreasing of the rel- 
ative weight of vertex corrections when approaching the 
split-band limit. Moreover, it displays a spurious peak 
being a consequence of using the local irreducible CPA 
vertex A in the expansion around the mean-field solution. 
When rewriting the expansion in terms of the full local 
CPA vertex 7 this peak vanishes. This peak is a rem- 
iniscence of a singularity in the irreducible vertex \ RR 
for A s» 1.5883 that is, however, compensated in the full 
vertex j RR , This fact gives evidence that one has to use 
the full local vertex in the expansion around mean-field 
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FIG. 4: (Color online) Ratio of the modulus of the vertex 
correction and the mean-field conductivity for conducti vity 
from Eq. (J28J), Acta, from Ref. 16, Aa^,, and from Ref. [lj, 
Acta/f. See the discussion in text. 



FIG. 5: (Color online) Vertex corrections to the conductivity 
calculated from Eq. I)28[). Acta, first-order approximation in 
7 to this expression, Acto, and the leading-order correction 
from the high-dimesional vertex from Eq. (|16[) . Acti. 



using the off-diagonal one-electron propagator as a small 
parameter. Otherwise we have no guarantee that the 
singularity from is completely compensated. 

A simplification in the momentum dependence of the 

integral kernel A from Rq. (|2Tj) is part of the non- 
perturbative approximation for the conductivity with 
vertex corrections, Eq. (|28|) . We replaced a product of 
two off-diagonal one-electron propagators by its average 
over the Brillouin zone. In Fig. [5] we compared the ver- 
tex correction from Eq. ()28[) with its second order ex- 
pansion in 7 (first order vanishes), Aoo, and the exact 
first-order vertex correction from Eq. ((4]) , Acti . The lat- 
ter two do not differ much in the strong-disorder limit, 
near the split-band transition where the vertex correc- 
tion is comparable with the Drude term. For weak dis- 
order, where the influence of the pole in function S(q) is 
smaller, approximate non-perturbative and perturbative 
vertex corrections almost equal and differ from the exact 
one. The impact of the vertex correction on the total 
conductivity is, however, negligible there. 

The primary objective of the presented construction 
was to rewrite and approximate the Kubo formula for 
the conductivity in such a way that two-particle vertex 
corrections to the one-electron Drude contribution never 
turn the total conductivity negative. Since the approxi- 
mate construction uses the full two-particle vertex, there 
is no obstacle to use this scheme as a global approxi- 
mation for one and two-particle quantities. One such 
a function of importance is the electron-hole correlation 
function $> RA . It is singular in the low-energy limit if 
the Ward identity is obeyed on a finite frequency inter- 
val around zero. Then also an Einstein relation couples 
the zero-temperature conductivity with the diffusion con- 
stant, a = e 2 pFD- As we argued, validity of the full scale 
Ward identity cannot be reached in theories distinguish- 
ing electron-electron and electron-hole scatterings. We 
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FIG. 6: (Color online) Scaling factor 
weight of the diffusion pole A -1 . 



from Eq. ([34]) and the 



can only guarantee validity of the Ward identity in the 
static limit lj — > 0. We introduced a reduction factor cf> 
with which we secured the existence of the diffusion pole 
in the electron-hole correlation function without the ne- 
cessity to correct the self-energy. We plotted in Fig.[5]this 
factor together with the weight of diffusive (extended) 
states A~ l . The reduction factor (j> does not differ much 
from unity even in the strong-disorder limit. The cor- 
rections to the self-energy induced by non-locality of the 
irreducible electron-hole vertex that had to be included 
in a fully self-consistent approach are hence not of prin- 
cipal relevance. As expected, the weight of the extended 
states vanishes at the split-band limit, A -1 = 0. 

The approximate conductivity is compared with those 
calculated via the Einstein relation from the diffusion 
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FIG. 7: (Color online) Conductivity and bare (D°) and renor- 
malized (D — D° /A) diffusion constants multiplied by density 
of states at the Fermi energy pF- 

constants D° and D = D°/A is plotted in Fig. [7] We 
can see that in the strong-disorder limit the conductiv- 
ity from the renormalized diffusion constant D matches 
the approximated one from the Kubo formula, while the 
conductivity from the bare diffusion constant D° is quite 
off for all disorder strengths. Note that D° > at the 
split-band limit. Conductivity a a and PfD° /A differ for 
intermediate disorder strengths but not (qualitatively) 
significantly. In the weak-disorder limit all three quanti- 
ties share the same asymptotics. 

VI. CONCLUSIONS 

We presented in this paper an approximate non- 
perturbative scheme for including non-local vertex cor- 
rections to the mean-field vertex from the coherent poten- 
tial approximation. The construction is a partial resum- 
mation of an expansion in the off-diagonal one-electron 
propagator from CPA that includes the leading-order 
non-local correction to the irreducible electron-hole irre- 
ducible vertex. After a simplification in the momentum 
dependence of the integral kernel of the Bethe-Salpeter 



equation with multiple scatterings we derived an explicit 
representation for the electrical conductivity, Eq. (|28|) . 
that remains non-negative irrespectively of how strong 
the disorder is. We explicitly demonstrated that this ex- 
pression produces reliable results for the electrical con- 
ductivity, in particular in the strong-disorder regime, and 
it qualitatively correctly mimics the behavior of the av- 
eraged conductivity in intermediate and weak-disorder 
regimes. An important aspect of the derived approach is 
that we made the approximation for the non-local part of 
the two-particle vertex compatible with the local mean- 
field one-electron self-energy. It means that it is not 
necessary to correct the self-energy by adding non-local 
terms in a cumbersome way to assess the averaged con- 
ductivity with vertex corrections. The only modification 
to be done in the electron-hole vertex is to rescale the 
local mean-field vertex 7 by a factor so that to restore 
the diffusion pole in the electron-hole correlation func- 
tion. This scaling effectively reduces the disorder effect 
in the regions where the vertex corrections to the conduc- 
tivity overweigh the Drude term and one needs to apply a 
non-perturbative approach to correct this deficiency. Our 
construction does exactly this with a minimum effort and 
no changes in the mean-field one-electron self-energy. It 
can hence be used to extend the standard mean-field one- 
electron approximation by including leading two-particle 
vertex corrections to the electrical conductivity without 
breaking consistency and the fundamental physical laws. 
The approximate construction is enough simple so that 
it may find application in the calculation of transport 
properties of real disordered materials beyond the model 
level presented here. 
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